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Abstract: Recently holographic techniques have been used to study the thermal properties 
of N = 2 super- Yang-Mills theory, with gauge group SU(N C ) and coupled to iV f <C N c 
flavours of fundamental matter, at large N c and large 't Hooft coupling. Here we consider 
the phase diagram as a function of temperature and baryon chemical potential fj, h . For fixed 
fj, h < N c M q there is a line of first order thermal phase transitions separating a region with 
vanishing baryon density and one with nonzero density. For fixed /i b > iV c M q there is no 
phase transition as a function of the temperature and the baryon density is always nonzero. 
We also compare the present results for the grand canonical ensemble with those for canonical 
ensemble in which the baryon density is held fixed [1]. 
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1. Introduction 

A large class of strongly coupled gauge theories can be studied using the gauge/gravity duality 
[2, 3]. In this context, the dynamics of a small number of flavours N f <C -/V c of fundamental 
matter in the gauge theory can be described by the dynamics of N f D-branes in the appropriate 
dual geometry At sufficiently high temperatures, for which the gauge theory is in a deconfined 
or plasma phase, 1 this geometry contains a black hole [4]. Since the condition N f <C A c ensures 
that the D-branes only perturb the background slightly, the dynamics of the fundamental 
matter in the deconfined phase corresponds to the dynamics of D-brane probes in a black 
hole background. 

The thermal properties of N f flavours of fundamental matter in SU(N C ) super- Yang- 
Mills (SYM) have been extensively studied using this holographic framework [5, 6, 7, 8, 9]. 
In these systems, the fundamental matter generically undergoes a first order phase transition 
at some critical temperature T iun . The low-temperature phase of the theory is described 
by 'Minkowski embeddings' of the probe branes (see fig. 1) in which the branes sit entirely 
outside the black hole [6, 7]. In this phase, the meson spectrum is discrete and exhibits a mass 
gap. Above the critical temperature T fun , the branes fall through the horizon and the gravity 
system is characterised by these 'black hole' embeddings. In this phase, the meson spectrum 
is continuous and gapless [7, 10, 11]. Thus, this large-iV c , strong coupling phase transition 
is associated with the melting or dissociation of mesons. As recently illustrated in related 
models [12], we emphasize that in theories that undergo a confinement /deconfinement phase 
transition at some temperature T doc < T fun , the phase transition for fundamental matter just 
described occurs as a separate phase transition. Instead, if T iun < T dcc , both phase transitions 
take place simultaneously. In the former case, mesonic states remain bound for a range of 
temperatures T dec < T < T fun despite the fact that the theory is in a deconfined phase. 

x The theory may or may not undergo a phase transition to a confined phase at low temperatures. 
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Figure 1: Possible D-brane embeddings in the black hole geometry. 

This physics is in qualitative agreement with that of QCD. Studies from lattice QCD [13] 
suggest that some mesons survive the deconfinement phase transition at T dec ~ 175 MeV and 
remain as relatively well-defined resonances up to temperatures of 2 — 3T dcc . It is therefore 
interesting to consider how this physics is modified in the presence of a chemical potential fj, b 
or nonzero density n b for baryon number. In the holographic framework, the introduction of 
a chemical potential [i b or nonzero density n b for baryon number corresponds to turning on 
the diagonal U(l) gauge field on the D-branes [1]. 

The phase diagram for the gauge theory at constant baryon number density n b was stud- 
ied in [1] . It was found that, for any nonzero value of the baryon number density n b , Minkowski 
embeddings are physically inconsistent. Simply stated, a nonzero density of baryons (or 
quarks) is dual to a worldvolume electric field, which in turn translates into a finite number 
of (fundamental) strings dissolved in the probe D-branes. Since the strings cannot simply ter- 
minate, it is not possible for the branes to close off smoothly above the horizon. Fortunately, 
for n b ^ there exist black hole embeddings where the D-branes intersect the horizon for 
all values of the temperature. This is in contrast with the case n b = 0, in which black hole 
embeddings do not exist below a certain temperature and Minkowski embeddings are required 
to describe the system [6, 7]. Focussing on the D3/D7 system, i.e., N f D7-brane probes in 
the black near-horizon geometry of N c D3-branes, it was shown in ref. [1] that the physics 
is essentially continuous around n b = 0. For small n b , the black hole embeddings mimic the 
behaviour of both the Minkowski and black hole branches of the n b = system, and a first 
order phase transition occurs between two different black hole embeddings - see ref. [1] for 
details. In a region with small n b , the (T,n h ) phase diagram exhibits a line of thermal phase 
transitions which terminates at a critical point (at finite n b ). However, it must be noted 
that the configurations below the phase transition were found to be unstable suggesting the 
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existence of a new, perhaps inhomogeneous, phase. Though the focus in ref. [1] was on the 
D3/D7 system, general arguments presented there suggest that the same physics occurs in 
other Dp/Dq systems, corresponding to gauge theories in different dimensions [14]. 

The above study of the thermodynamics with fixed n h was performed in the canonical 
ensemble. Of course, the thermodynamics can also be studied in the grand canonical ensemble, 
in which the chemical potential // b is kept fixed. This is the focus of this paper. As described 
above and clearly elucidated in [1], the restriction to black hole embeddings only applies when 
n b is fixed and nonvanishing. Hence, as we will describe, Minkowski embeddings play an 
essential role in describing the grand canonical ensemble. Again we consider only the D3/D7 
brane system in the following but expect similar results for other Dp/Dq brane systems [14]. 

Our main result is the phase diagram displayed in fig. 2, where 



(i.i) 



is the quark chemical potential and M tx M q is a mass scale defined in eqn. (2.14) below. As 
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Figure 2: Phase diagram: Quark chemical potential /x q /M q versus temperature T/M. The red line 
separates the phase of Minkowski embeddings (small temperatures, small fj, q /M q ) from black hole 
embeddings. Figure (b) zooms in on the region near the end of this line and also depicts the boundary 
of the region accessed by the black hole embeddings (green) and a small region (enclosed by the blue 
curve) where more than one black hole embedding is available for a given value of // q and T. The 
separation of the red and green curves would not be resolved on the scale of Figure (a). 



we show in section 2, for fixed /i q < M q black hole embeddings do not exist below a certain 
temperature. In contrast, for /U q > M q , black hole embeddings exist for all temperatures. On 
the other hand, Minkowski embeddings with a vanishing density of baryons, n b = 0, exist for 
all /i q up to a temperature T ~ T fun . Our study in section 3 of the free energy of the system as 
the temperature increases reveals that for /i q < M q there is a first order phase transition from 
a Minkowski embedding to a black hole embedding. In the field theory this is a transition 
characterised by the condensation of charge, namely by a jump from a phase with ra q = 
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to a phase with n q 7^ 0. The meson spectrum across this transition also changes: In the 
low-temperature Minkowski phase the meson spectrum is discrete and possesses a mass gap, 
whereas in the high-temperature black hole phase it is continuous and gapless. For larger 
chemical potential, /i q > M q , there is no phase transition as a function of the temperature. 

As described above, for // q < M q black hole embeddings do not exist below a certain 
temperature. The boundary of the portion of the phase diagram accessed by these embeddings 
essentially coincides with the line of phase transitions shown in fig. 2a. However, at higher 
resolution in fig. 2b, one can see these two boundaries deviate in a small region near p q ~ 0. 
This figure also illustrates a small region of the phase diagram where more than one black 
hole embedding is available for a given value of /x q and T. In particular, the latter contains 
the black hole embeddings that are thermodynamically unstable. As will be discussed below, 
this 'multi- valued' region of black hole embeddings was central to the analysis of the phase 
transition studied in [1] for the canonical ensemble with fixed n q . In section 4, we discuss the 
reconciliation of these previous results with those found here. 

While this paper was in the final stages of preparation, three related papers [15, 16, 17] 
appeared which have considerable overlap with our present investigation. 



2. Holographic framework 

Following [6, 7], the metric for N c black D3-branes in the decoupling limit may be written as 2 

(2.1) 



2 1 fuop\ 2 

ds = 2 VT) 
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- J -^dt 2 + fdx 2 
f 



I? 



+ [dp 2 + P 2 dn 2 ] , 



where p is a dimensionless coordinate and 



f(p) = l-^, f(p) = l + ^, L 4 =4irg s NJt 
This metric possesses a horizon at p = 1 with temperature 

m U 
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(2.2) 
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which is identified with the temperature of the dual J\f = 4 SYM theory. In turn, the 
coordinates {t, x} are identified with the coordinates of the gauge theory. The string coupling 
constant is related to the SYM 't Hooft coupling constant through 



A = gl M N c = 2vr 5s iV c . 



(2.4) 



In addition to the metric above, the D3-brane solution is characterised by a constant dilaton 
and a Ramond-Ramond field given by 

{uopff 2 



C0123 — — " 



4L 4 



(2.5) 



2 This metric is related to the standard presentation through the coordinate transformation (uop) = u 2 + 
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Introducing N ( D7-branes into the geometry above corresponds to coupling N ( fundamen- 
tal hypermultiplets to the original SYM theory [18]. Before the decoupling limit, the branes 
are oriented as 

0123456789 
D3: x x x x (2.6) 
D7: xxxxxxxx 

This configuration is supersymmetric at zero temperature, which ensures stability of the 
system. After the decoupling limit, the D7 branes wrap an S 3 of possibly varying radius 
inside the S 5 of the background geometry. Writing 

dfi| = d9 2 + sin 2 9dfl 2 3 + cos 2 9dxj> 2 (2.7) 

and defining \ = cos 9, the embedding of the D7-branes is described by 4> = 0, x = x(p)- 

The derivation of the equations of motion for the D7-branes' embedding x{p) an d the 
gauge field on their worldvolume At was discussed in [1] and we review a few salient details 
here. The (Dirac-Born-Infeld) action of the D7-branes is 



j D7 = -iv f r D7 J dt d 3 xd P dn 3 ^^ff(i - X 2 )^i- X 2 + P 2 X 2 _ _ X 2 } i t 2 ; (2 8) 

where the dot denotes differentiation with respect to p and we have introduced a dimensionless 
gauge field 

A t = ^A t . (2.9) 

Note that At corresponds to a radial electric field. The corresponding equation of motion 
(eqn. (2.11) in [1]) dictates that asymptotically 

* = *-2^ + -' (2 - 10) 

where the constant /U q is the quark chemical potential and d is a dimensionless constant related 
to the vacuum expectation value of the quark number density operator through 

n q = ^N f N c V\T 3 d. (2.11) 

The equation of motion for \ (eqn. (2.17) in [1]) implies that the branes' profile behaves 
asymptotically as 

m c .„ 

X = - + - + •••> 2 - 12 

p p A 

where the dimensionless constants m and c are proportional to the quark mass and condensate, 
respectively [6, 7]. In particular 

M 

m=— , (2.13) 
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where the scale 

M = ^ = ^ (2.14) 
2vr K ! 

is (up to a factor) the meson mass gap M gap in the D3/D7 brane theory at zero temperature 
[20]. 

Let us now consider what types of embeddings are possible if // q / 0. As shown in [6, 7], 
Minkowski embeddings only exist for T/M < 0.77, i.e., for sufficiently low temperatures or 
sufficiently large quark masses. Moreover, these embeddings are physically consistent only if 
n q = [1]. However, the chemical potential can take any value, since a Minkowski embedding 
with At = /i q is a solution of the equations of motion if the same embedding is a solution 
with At = 0. Thus Minkowski embeddings describe a phase with a possibly nonzero chemical 
potential but vanishing charge density. We therefore expect (and we will confirm below) that 
these embeddings are thermodynamically preferred if /i q < M c , where M c is the constituent 
or thermal quark mass at a given temperature. As shown in fig. 1, this mass is proportional 
to the distance between the tip of the D7-branes and the horizon, in contrast to the bare 
quark mass M q , which is proportional to the asymptotic distance between the D7-branes and 
the D3-branes [7, 19]. 

If n q / then only black hole embeddings are physically consistent, and these exist 
for all values of T/M. Regularity requires the gauge field to vanish at the horizon, so the 
chemical potential can be expressed in terms of d, T and M q by integrating Gauss' law from 
the horizon to the boundary as [1] 



/(1- X 2) ,6/3(1 _ x 2)3 +8d 2 



A key observation now is the fact that, for sufficiently small T/M, n q approaches a nonzero 
value in the limit n q — > 0. This is illustrated in fig. 3, which displays /x q /M q as a function 
of T/M for decreasing values of d (from top curve to bottom curve) all the way down to 
d = 10 -4 . Taking d to be even smaller than this value does not yield values of jtx q appreciably 
smaller on the scale shown in the figure. In other words, black hole embeddings only cover 
the region above the bottom curve in fig. 3. This observation is important because it means 
that Minkowski embeddings are necessary to describe the region below this curve. 

The curves in fig. 3 were obtained by numerical integration of eq. (2.15). However, 
analytic results can be obtained in the limiting cases of large and small T/M. In the limit 
T/M — ► with fixed d, it was shown in [1] that /U q — ► M q regardless of the value of d. This 
is reflected in fig. 3 by the fact that all curves meet on the vertical axis at /U q /M q = 1. In 
the opposite limit, M/T — > 0, the embedding becomes simply \{p) = 0- Substituting this 
into (2.15) one obtains an expression that can be easily integrated. In particular, one easily 
verifies in this case that 
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Figure 3: Chemical potential /z q /M q versus T/M for various values of d, increasing from bottom up: 
d = 1CT 4 , 1CT 2 , 1, 10, 104, 4298, 380315. The J= 10~ 4 , 10~ 2 lines are virtually coincident. 

for small d, and consequently /i q —* as d —* 0. 

To summarise, we conclude that the region of small // q /M q and small T/M (the trian- 
gular region below the bottom line in the lower left corner of fig. 3) must be described by 
Minkowski embeddings. Similarly, the region with T/M > 0.77 must be described by black 
hole embeddings. In order to decide what the favoured embedding is in the rest of the pa- 
rameter space, i.e., in the region in which both Minkowski and black hole embeddings exist, 
we must determine which embedding minimises the free energy of the system. 

3. Thermodynamics 

From the previous section it is clear that both Minkowski and black hole embeddings play a 
role in the grand canonical ensemble. As both embeddings exist for certain values of temper- 
ature and chemical potential, it is necessary to study the thermodynamics to determine the 
physically preferred configuration. To this end, we use the standard technique of Wick rotat- 
ing [21] the time direction. Then, the leading contribution to the free energy is determined 
by evaluating the Euclidean action. As we are interested in the behaviour of the fundamental 
matter, we only consider the action for the D7-branes. 

Of course, the Euclidean action for the D7-branes is given by Wick-rotating the action in 
(2.8). Evaluating this action leads to a formally divergent result. However, the divergences are 
removed by introducing a finite-radius ultraviolet cut-off and a set of boundary counterterms 
to renormalise the action [23]. This holographic renormalisation was discussed in detail in 
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[1, 6, 7] and the result, from ref. [1], is 



& = G < m > 



[G>L 



m 



2\2 



4mc] , 



where G{m) is the convergent integral 



G(m) 



dp P 3 ff(l 




X 2 + P 2 X 2 



2 /, 



1 - X 2)i| 



Q 2 

p + m p 



N is a normalisation constant given by 



= —\N C N C T 3 , 
32 ' 



(3.1) 



(3.2) 



(3.3) 



and /9 min is the minimum value of p for the embedding of interest. 

This Euclidean action of the D7-branes is identified with their Gibbs free energy W(T, /j, q ) 
via W = TI E [1]. Of course, the thermodynamically preferred D7-brane embedding is 
that which minimises the free energy and hence we evaluated the free energy as a func- 
tion of temperature for fixed p q /M q . Representative plots of the free energy (scaled by 
J\f = XN f N c M 4 /32) versus temperature are provided in figs. 4 and 5. The results for 
p q /M q = 0.0108 (fig. 4) are typical for small p q /M q , displaying the classic 'swallow tail' 
shape and closely resembling the /i q = results (shown in fig. 5 of [7]). For larger fijM v the 
'swallow tail' shape disappears - see fig. 5. For p q /M q < 1, the two branches still cross at 
some temperature and this, of course, determines the temperature of the phase transition - 
see fig. 5a. As p q /M q increases towards 1, the temperature of the phase transition decreases. 
The results in fig. 5b are typical for large p q /M q : The two branches never cross and the free 
energy is always lower on the black hole branch. 




0.760 0.762 0.764 0.766 0.768 U.770 0772 



Figure 4: Free energy versus temperature for /i q /M q = 0.0108187. The blue dotted line (red solid) 
represents the Minkowski (black hole) branch. The vertical line marks the temperature of the phase 
transition. 

By varying fj, q /M q and examining the free energy, the phase diagram in fig. 2 was mapped 
out. For [i q jM q < 1, there is a phase transition corresponding to a discontinous jump from 
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(a) njM q = 0.8 (b) /x q /M q = 2 

Figure 5: Free energy versus temperature for (a) /x q /M q = 0.8 and (b) ^t q /M q = 2. In (a), the vertical 
line marks the temperature of the phase transition. 

a Minkowski to a black hole embedding. As /i q /M q increases towards one, the temperature 
of the phase transition decreases. For larger chemical potentials, // q /M q > 1, there is no 
phase transition, with black hole embeddings thermodynamically favoured for all temper- 
atures. Other thermodynamic quantities can be computed from the free energy, e.g., the 
entropy S = (dW/dT) . Discontinuities in thermodynamic quantities (e.g., quark conden- 
sate, entropy density) at the phase transition indicate that the phase transition is first order. 
As illustrated in fig. 5, however, the size of the discontinuity is becoming smaller and smaller 
as we approach ^ q /M q = 1. In fact, recent analytic work [17] indicates that this is a critical 
point where the phase transition becomes second order. For [i q /M q < 1, the low temperature 
phase is described by Minkowski embeddings and hence n q = below the phase transition. 
After the phase transition, n q > and it increases as the temperature increases. Similarly, 
for /i q /M q > 1, n q increases as the temperature increases. 

Finally, we turn briefly to the thermodynamic stability of the system. The conditions 
for stability can be stated in different ways. Generally they can be phrased as the restriction 
that the Hessian of the state function, here the Gibbs free energy, is positive semi-definite. A 
simpler approach is to examine the necessary (but not sufficient) conditions: 



dS_ 
df 



<9/x q 



>0, (3.4) 

T 



which ensure that the grand canonical ensemble is stable against fluctuations purely in energy 
or charge. In ref. [1], it was stated that in the canonical ensemble the first requirement with 
fixed n q (rather than fixed /it q ) appears to be satisfied everywhere, however, more recently 
ref. [15] reported that there are some regimes where the latter condition fails to hold. At 
present, we have found no indications that the grand canonical ensemble suffers such an 
energetic instability. Further, ref. [1] found that the second condition does not hold for some 
black hole embeddings. 
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These unstable configurations can be found by examining fig. 6, which illustrates as a 
function of T and n q for the black hole embeddings - note that d/m 3 oc n q /M q is independent 
of T. The figure only covers a small region near the phase transition found in the canonical 
ensemble. In this regime, the surface /z q (T, n q ) displays an interesting 'fold' and in the phase 
transition described in [1], the system jumps between configurations on the top and bottom 
of this fold. One can see that the unstable configurations discussed above correspond to the 
top of the fold illustrated in fig. 6. These unstable embeddings then play a central role in the 
phase transition found for the canonical ensemble. In fig. 2, these unstable configurations 




T/M 

Figure 6: Three-dimensional plot of the chemical potential, the temperature and the charge density 
determined by black hole embeddings. 

appear near the top of the multi-valued region enclosed by the blue curve and hence they 
could in principle play a role in this regime for the grand canonical ensemble as well. However, 
it turns out that minimizing the Gibbs free energy always picks out either a stable black hole 
or Minkowski embedding in this regime. Hence the system never suffers from this electrical 
instability. This result poses an apparent paradox though since (in the infinite volume limit, as 
applies here) the canonical and grand canonical ensembles should be equivalent [22], whereas 
here they seem to display different behaviours. We address this issue in the discussion below. 
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4. Discussion 

A universal, first order, thermal phase transition in Dp/Dq systems was identified in [6], 
characterised by a jump of the Dq-branes from a Minkowski to a black hole embedding in the 
Dp-brane geometry. In the dual gauge theory, the transition is associated with the melting 
or dissociation of mesons. The analysis of ref. [6] was extended to include a nonzero quark 
density in ref. [1] , where it was found that Minkowski embeddings are physically inconsistent 
for any n q ^ 0. Despite this, a first order thermal phase transition, in this case between two 
black hole embeddings, was found to occur for densities below a certain critical density. 

In this paper we have examined the thermal behaviour of the D3/D7 system as a function 
of the quark chemical potential /x q , as opposed to the quark density n q . As described in the 
introduction and clearly elucidated in [1], the restriction to black hole embeddings only applies 
when n q is fixed and nonvanishing 3 and so no longer applies in describing the grand canonical 
ensemble. In particular then, we have found that, for /i q < M q , the low-temperature phase 
is described by Minkowski embeddings with n q = 0. If the temperature is kept fixed and the 
chemical potential is increased to a critical value, or viceversa, the system undergoes a first 
order phase transition characterised by a jump to a black hole embedding with n q ^ 0. Thus 
the phase transition is characterised by the condensation of charge. The meson spectrum 
for the low-temperature phase was calculated in ref. [7] and is discrete and exhibits a mass 
gap. In the high-temperature phase there are no stable mesons and we expect a discrete 
spectrum of quasinormal modes [10, 11]. Hence, as with the (i q = case [6, 7], a dramatic 
feature characterising the phase transition is the melting of mesons. Unlike in the fi q = case, 
however, there is a finite density of quarks in the high-temperature phase, and this density 
increases with temperature. Fig. 2a shows that the temperature of the phase transition 
decreases as fi q /M q increases towards one. 

For /x q > M q , there is no phase transition: In the gravity system, there are black hole 
embeddings of the D7-branes corresponding to all temperatures. Computing the free energy 
of these configurations and comparing to the results for the Minkowski embeddings reveals 
that the black hole embeddings are thermodynamically preferred. It would be interesting to 
examine the effect of a nonzero chemical potential on the transport properties of the D3/D7 
plasma in this regime by extending the analysis of refs. [8, 10, 24, 25]. 

In this black hole phase we expect the meson spectrum to be continuous and gapless. 
The true quasi-particle excitations of the system, if any, correspond to 'quasinormal modes' 
of the D7-brane worldvolume fields [11]. In a regime in which the imaginary part of the 
quasinormal frequency is small, these excitations can be interpreted as quasi-particles in 
the dual quark-gluon plasma. Such long-lived quasi-particle excitations yield particularly 
striking resonances in the corresponding spectral functions. Such resonances were found for 
the D3/D7-brane system with n q = in [10, 25] but these were most pronounced for unstable 
black hole embeddings. Hence, it would be particularly interesting to study the spectrum of 

3 Refs. [15, 26] attempt to construct Minkowski embeddings with n q 7^ but we believe these are unphysical 
for the reasons described in [1] . 
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quasinormal modes [10, 11] or spectral functions [10, 25] in detail in the present context since, 
as described in [1], we expect to find a class of almost stable meson excitations in the regime 
of low temperature and low quark density where the black hole embeddings should be stable. 

Let us expand on this point: In this regime of low temperature and low quark density, the 
D7-brane connects to the event horizon with a long narrow spike. This black hole embedding 
then mimics the behaviour of a Minkowski embedding with (a gas of) fundamental strings 
stretching down to the horizon [1]. Hence the quasinormal spectrum should include long-lived 
modes that mimic the spectrum of (stable) mesons on the corresponding Minkowski embed- 
ding. It is relevant here to recall that the Minkowski phase also contains stable unconfined 
quarks described by strings stretched between the brane and the horizon. We note that these 
quarks also have unstable excited states described by the quasinormal modes on the strings. 
Presumably the quasinormal spectrum of the black hole embedding must also include modes 
which mimic this behaviour. On the black hole embedding, these modes describing 'quark 
excitations' would have support primarily on the spike connected to the horizon, in contrast 
to the 'mesonic' modes whose support would be primarily above the spike. Given the stan- 
dard intuition for holographic energy scales [27], one may expect T < Re(a->) < M q /\/A for 
the 'quark excitations'. In contrast, the imaginary part of the longest-lived modes is of order 
Im(w) ~ \ r \T 2 /M q < T in the regime considered here [19]. Hence it is possible that these 
modes may represent an independent class of quasiparticle excitations in the quark-gluon 
plasma. 

In the Minkowski phase the spectrum contains mesons and quarks, and both sets of 
excitations are absolutely stable in the large- N c , large-A approximation considered here. In 
view of the presence in the spectrum of stable, quark-number-carrying states, it may seem 
surprising that the physics in the Minkowski phase is completely independent of /x q , in par- 
ticular that the quark density vanishes in this phase. After all, at finite temperature one 
expects a gas of quarks and anti-quarks to be present, and thus a finite chemical potential to 
induce a charge imbalance. The resolution of this puzzle comes from the fact that the quark 
mass is of order M q , whereas the maximum temperature for which Minkowski embeddings 
are thermodynamically preferred is set by the meson mass and is therefore parametrically 
smaller: T ~ T fun ~ M q /vA. Thus the abundance of quarks and anti-quarks in the thermal 
gas is proportional to the Boltzman factor 

e -M q /r^ e -VA_ 

This is exponentially suppressed at strong coupling, and hence is presumably invisible on the 
gravity side at leading order in the large-A, large- N c approximation under consideration. 

One may also wonder whether baryons may play a role in the Minkowski phase. From 
the gauge theory viewpoint, we do not expect this to be the case, since the baryon mass scales 
as N c and hence the baryon abundance at finite temperature is suppressed by exp (—N c ). On 
the gravity side, at strictly zero temperature, stable baryons are represented by D5-branes 
wrapping the S 5 part of the geometry [28] and connected to the D7-branes to induce an 
electric flux corresponding to N c fundamental strings [20]. However, our intuition is that 
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upon introducing a nonzero temperature, the gravitational attraction of the black hole causes 
the D5-branes to collapse and they are pulled behind the event horizon leaving behind just 
a black hole embedding of the D7-branes. We have performed some numerical investigations 
that seem to support this conclusion, but our results remain preliminary. Specifically, we 
have investigated the possible existence of tear-drop configurations, 4 similar to those of [29]. 5 
These configurations would end on a cusp above the horizon that could be connected to 
a cuspy D7-brane configuration along a number of fundamental strings. The tension of the 
strings would then overcome the tension of the branes [1] and presumably force the D5-branes 
and the D7-branes to connect to each other, resulting in a smooth configuration. However, 
our preliminary results indicate that there are no static configurations of this type. It is 
interesting to imagine a dual description of the dynamical collapse of the suggested tear-drop 
configurations. One would begin with a baryon as an excited state of bound quarks which 
are then dispersed as free quarks by the thermal bath of the strongly coupled gauge theory. 
In any event, we re-iterate that our numerical investigations here are only preliminary. 

In the limit T — > or M q — ► oo, the phase transition from a Minkowski to a black hole 
embedding happens exactly at /x q = M q . One can see this as arising from the fact that the 
D3/D7 system is super symmetric at T = 0. In this limit the force exerted on the D7-branes 
by the background sourced by the D3-branes vanishes and the D7-branes are not deformed 
[20]. This implies that the constituent quark mass exactly coincides with the bare quark 
mass, M c = M q . It is therefore not surprising that the phase transition takes place when the 
chemical potential reaches this mass. The constituent quark mass for T > (on Minkowski 
embeddings with n q = 0) was calculated in [7, 19]. Fig. 7 shows the ratio M c /M q (blue, 
dotted curve) as a function of temperature, as well as the value of the ratio /u q /M q (red, 
continuous curve) at which the phase transition takes place. It is remarkable that the two 
curves are almost coincident except for relatively small values of /x q , i.e., towards the higher 
temperature end of the line of phase transitions. A similar result was found recently in [16] 
for a similar holographic gauge theory. We note that the two curves in fig. 7b begin to diverge 
significantly at T ~ .758M, which coincides with the temperature where the line of phase 
transitions enters the multi- valued region in fig. 2b. This agreement seems better than one is 
in principle entitled to expect. That is, fig. 7 illustrates that the chemical potential required 
to introduce quarks into the gluon plasma essentially matches the energy of a 'free' quark. 
We have put 'free' in quotes because the constituent quark mass M c does take into account 
the interactions of an individual quark with the strongly coupled plasma of adjoint fields. 
However, the agreement in fig. 7 then indicates that the interactions between the quarks 
themselves are largely negligible at the phase transition. 

4 One might also consider cone- like configurations in which the D5-brane has already partially fallen through 
the event horizon. These would describe a mixed phase of baryons and unbound quarks, if they were stable, 
but our intuition is that they are also unstable. 

5 Note that the D5-brane configurations studied in [29] do not correspond to baryons in the four-dimensional 
M = 4 SYM theory at finite temperature. The gravitational background in [29] corresponds to the AdS soliton 
rather than an AdS black hole, so the baryons considered in that reference are baryons in a confining phase of 
the gauge theory produced by compactifying one of the dimensions. 
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Figure 7: Comparison of the ratio M c /M q (blue, dotted curve) with the ratio /x q /M q (red, continuous 
curve) at which the phase transition from a Minkowski to a black hole embedding takes place. The 
two curves essentially coincide on the scale of the Figure (a). 



This result may further suggest that the density of quarks should be small just above 
the phase transition. Examining fig. 8, we see that this statement seems to be correct in that 
the dimensionless ratio r = d/m s is indeed much less than one. Further we note that r is 
largest towards the higher temperature end of the line of phase transitions, i.e., roughly in 
the same region where the two curves in fig. 7b deviate from one another. However, using 
eqs. (2.11-2.14), we find 

(M \ 3 

n q = V2r\N t N c n clit with n crit = (^J . (4.2) 

Above we have defined n crit keeping in mind that studies of structure functions [30], as well as 
mapping of a quark's disturbance of the adjoint fields [31], both indicate that the size of (the 
glue cloud around) an individual quark is roug hly {MjVX)- 1 (at T = 0). Hence we might 
expect that the quarks would begin to interact significantly when the quark density reaches 
n crit . Therefore, even though r may be of the order .0001 (as at the phase transition in fig. 8), 
eq. (4.2) indicates that n q 3> n crit because we are in the limit of large A and large N c . Hence 
the good agreement between ^ q and M c suggests the surprising result that the quark-quark 
interactions have a negligible effect even at parametrically large quark densities. 

Let us comment on the supergravity description of these results. The constituent quark 
mass is defined as the energy of a test fundamental string stretching between the tip of the 
D7-branes and the horizon [7, 19]. In the absence of supersymmetry (i.e., at T > 0), this 
energy need not coincide with the energy of a 'spike' emanating from the D7-branes down 
to the horizon. This spike in the D7-brane geometry extends uniformly in all of the gauge 
theory directions. Hence this spike accounts for the effects of the interactions of the quarks 
with both adjoint fields and other quarks. In other words, this geometry, which represents 
the presence of a macroscopic string density, has an energy which receives contributions not 
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accounted for by the test string calculation, such as (for example) the deformation of the 
D7-branes. 



0.0035 




T/M 

T/M 

(a) (b) 

Figure 8: Charge density at the phase transition in the grand canonical ensemble (red). Figure (b) 
shows in blue the line of phase transitions identified in the canonical ensemble at fixed n q [1]. Further 
the shaded region enclosed by the green and blue curves corresponds to the unstable region identified 
in [1], in which (9^ q /9n q )T < 0. 

Before closing we would like to comment on the relation between the present results and 
those presented in [1]. There for the canonical ensemble with fixed n q , we found a first order 
phase transition that only existed for small charge densities and temperatures near T fun , the 
critical temperature identified at n q = 0. Here for the grand canonical ensemble with fixed 
fiq, we found a first order phase transition extending from (/i q ,T) = (0,T fun ) to (M q ,0), as 
illustrated in fig. 2a. The difference was further highlighted at the end of section 3 with the 
observation that certain unstable black hole embeddings play a role in describing the phase 
diagram for the canonical ensemble - the corresponding region is illustrated in fig. 8b - but 
they are not physically relevant for the present analysis of the grand canonical ensemble. 
These discrepancies pose a paradox since (in the infinite volume limit under consideration 
here) the canonical and grand canonical ensembles should be equivalent [22]. 

Before resolving the inconsistency, let us elaborate on what the equivalence stated above 
entails at a practical level. In the grand canonical ensemble we fix T and /i q . Minimizing 
the Gibbs free energy then determines the charge density n q = n q (T, /z q ). In the canonical 
ensemble, we fix T and n q and minimize the Helmholtz free energy. This determines the 
chemical potential /i q = /^ q (T, n q ). However, both of these surfaces, (T, /i q , n q (T, // q )) grand and 
(T, Hq(T, n q ), n q ) canon , should coincide. 6 That is, a given T and // q in the grand canonical 
ensemble determine a particular value of n q . Then with this value of n q and the same 

6 While this statement is correct for practical purposes, there are, of course, technical subtleties, e.g., at 
the line of phase transitions where dfi q /dn q = 0. Similar issues arise at n q = 0, where 9n q /cfyi q = in the 
low-temperature phase with the present large-A, large- iV c approximation. 
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temperature T, thermal equilibrium in the canonical ensemble will yield the same value of /Lt q 
as before. 

Now the resolution of our inconsistency can be found by examining fig. 8 which illustrates 
the quark density on the black hole embedding along the line of phase transitions in fig. 2. 
Recall that on the other side of the phase transition, the Minkowski embeddings have n q = 
(or d = 0). Hence if we map each configuration in the phase plane in fig. 2a to the (n q , T) 
plane as in fig. 8, we find a gap between the n q = axis and the red curve illustrated 
in the second figure. Hence we are left to explain how the system accesses these values 
of n q and T. Here we can use the intuition that we might derive from, say, the liquid to 
gas transition of ordinary water. That is, we expect that this intermediate regime is filled 
in by an inhomogeneous mixture where the two phases (with the same chemical potential) 
coexist. Now the key point, shown in fig. 8b, is that the phase transition identified for the 
canonical ensemble in [1], as well as the region where unstable embeddings were favoured in 
this ensemble, lies entirely within this intermediate regime where we have just argued that 
the system should be inhomogeneous. Since the analysis of [1] was restricted to only consider 
the homogeneous configurations, it misidentified the correct ground state of the system in the 
small region around the phase transition. 

Of course this is not a surprise since, as discussed in [1], the electrical instability of the 
configurations below the phase transition there already pointed towards an inhomogeneous 
phase playing a role in the correct phase diagram. The present discussion gives a clearer 
picture of the nature of this inhomogenous phase and where it is relevant. Let us add that 
this clarification also highlights a shortcoming of the holographic explorations of the phase 
diagram of strongly coupled gauge theories here and, for example, in [1, 15, 16, 17, 26, 33]. 
That is, they rely on constructing simple and highly symmetric configurations of D-branes 
and supergravity fields - in particular, configurations describing homogeneous phases of the 
gauge theory. These symmetries are imposed by hand and, since it is difficult to account 
for deviations towards less symmetric configurations, the latter are typically ignored. As 
highlighted by the present discussion, this simplistic approach may fail even in relatively 
simple circumstances, and so one must be cautious in extracting conclusions based on highly 
symmetric configurations. 

Let us finish with some comments on the possible comparison of the phase diagram 
obtained here with that of QCD. In addition to the usual caveats at zero density, there are 
two additional caveats that are particularly relevant at finite chemical potential. The first one 
is the fact that baryon number in QCD is only carried by fermionic fields, whereas in the type 
of systems we have studied it is also carried by scalar fields in the fundamental representation. 
This feature is absent in the D4/D8/D8 system of Sakai and Sugimoto [32], whose spectrum in 
the fundamental representation consists only of fermions. Interesting analyses of this system 
at finite baryon number chemical potential have recently appeared [33], and further studies 
will likely be a fruitful avenue for the future. The second caveat is the fact that the very 
existence of many of the expected phases of QCD at finite chemical potential depend crucially 
on the fact that N c = 3 (see e.g., [34, 35]). It is therefore unclear to what extent the large- 
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N c limit may provide a useful approximation at finite chemical potential; some optimistic 
arguments, however, can be found in [36]. In any case, given the lack of other analytical 
methods and the difficulties of simulating QCD at finite chemical potential on the lattice, 
we think it is interesting to explore possible insights coming from the holographic duals of 
large- N c gauge theories. 
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